BCS-BEC crossover with a finite-range interaction 
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We study the crossover from BEC to BCS pairing for dilute systems but with a realistic finite- 
range interaction. We exhibit the changes in the excitation spectrum that provide a clean qualitative 
distinction between the two limits. We also study how the dilute system converges to the results 
from a zero-range pseudo-potential derived by Leggett. 
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I. INTRODUCTION 

Ultracold atomic gases provide an experimental play- 
ground for testing pairing phenomena, due to the 
ability to control the inter-atomic interactions via a 
magnetically-tuned Feshbach resonance^. Over the past 
year, tremendous progress has been made in realizing 
the crossover from the Bose-Einstein condensation (BEC) 
of diatomic molecules to the Bardccn-Coopcr-Schricffcr 
(BCS) limit of weakly-bound Cooper pairs in ultracold 
gases of fermionic atoms^i But there is much debate 
over where the crossover between the BCS and BEC 
regimes occurs, following recent experimental studies on 
fermionic condensates^. Therefore, it is timely to per- 
form a thorough theoretical investigation of what the sig- 
nature of the crossover is and how it might be observed 
experimentally, using a realistic potential with a finite 
range. 

Much of the theoretical work on systems of two types 
of fermions interacting via an adjustable, attractive po- 
tential have focussed on interactions that are governed 
by a single parameter, namely the s-wave scattering 
length ay^. Such a description is valid provided we 
have | do | 3> (r) and kp(r) -C 1, where (r) is the width 
of the potential and kF is the Fermi momentum of the 
non-interacting system, so that the only independent di- 
mcnsionless variable in the problem is a^kp- Thus, it 
is only applicable to dilute systems like the ultracold 
atomic gases, not the high density situation found in 
conventional superconductors^. The potential supports a 
2-body bound state for (cio^f) -1 > 0, but this molecu- 
lar state passes through zero energy and vanishes into the 
continuum at (aofcp) -1 = 0, the position of the Feshbach 
resonance. As such, the BCS and BEC limits correspond 
to (aokp)^ 1 — > — oo and (dofcp ) - * +oo, respectively. 

The distinction between the BEC and BCS regimes can 
be made sharper—: the qualitative boundary between the 
BCS and BEC ground states occurs when the chemical 
potential /x reaches 0, marking the disappearance of a 
Fermi surface, and this coincides with the appearance 
of a bound state at (oo/cf) 1 = only when the den- 
sity is zero. At the point where there is no longer a 
defined Fermi surface, there are qualitative changes in 



the quasi-particlc excitation spectrum: the momentum 
corresponding to the minimum in the excitation spec- 
trum shifts from finite momentum in the BCS limit to 
zero momentum in the BEC limit. At the same point, 
the value of the excitation gap goes from A in the BCS 
limit to -\//i 2 + A 2 in the BEC limit, where A is the gap 
parameter. We shall use these changes in the excitation 
spectrum as the "definition" of the crossover point. 

There have also been studies of the BCS-BEC crossover 
using Coulomb potentials in the context of excitonst 
Here, the interaction is kept fixed and the density varied 
to achieve a smooth transition from the high density BCS 
limit to the dilute BEC limit. Note that it is not possible 
to obtain a density-driven crossover using an interaction 
that is only characterized by the scattering length, since 
this picture breaks down at high densities. 

In this paper, we shall study the BCS-BEC crossover 
at zero temperature using a realistic Gaussian potential 
and compare with results from single-parameter poten- 
tials. We use a mean-field variational wave function to 
describe the crossover, which is expected to give an accu- 
rate description of the low temperature behavior, though 
not, of course, the finite temperature transition. 

The paper is organized as follows: We describe the 
model and variational scheme in Sec. [H] We discuss the 
characteristics of the wave functions and excitation spec- 
tra in Sec. IIIII In Sec. IIVI we compare our results with 
Leggett's zero-range pseudo-potential result which is usu- 
ally used to describe the mean-field limit for dilute Fermi 
systems. We conclude in Sec. [V] 



II. FORMALISM 

We determine the ground-state wave function | $) using 
conventional mean-field methods, briefly reviewed here. 
We consider the Hamiltonian 

H = Yl £k a ^ akl + J2 V i a kT a pi a P-qi ak +iT ' C 1 ) 

k,i k.p.q 

where k and i denote momentum variables and spin 
states {T)j}j respectively, and ek = ft 2 k 2 /2m. Since ex- 
periments on atomic gases are carried out at low energies, 
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we will restrict ourselves to the simple case of an s-wave 
pairing interaction. 

We calculate the ground state by minimizing the free 
energy 



F = ($\H-(j,N\$) , 



(2) 



where N = i aj^aici is the total number operator and 
fi is the chemical potential. The standard BCS varia- 
tional wave function is given by 



|$) = M exp (^2 



-4^-14 1°)' ( 3 ) 



and it smoothly interpolates from the BCS to BEC lim- 
its, giving an accurate description of the crossover. Here, 
the BCS parameters u k , «k only depend on k = |k| in the 
s-wave approximation, and Af = J} k u k is the normaliza- 
tion constant, such that ($|$) = 1. With this ansatz, 
the free energy becomes 

F = 2 Ok - A*) pk + V v - k K* kKp , (4) 

k k,p 

where the normal and anomalous densities are 



Pk = ($|a kT a kT |$) = kk| 2 , 



(5) 



and 



K k = ($|a_kiakTl ( S) = vlu u , (6) 

respectively, subject to the normalization condition, 

\u k \ 2 + M 2 = 1. (7) 

Without loss of generality we can choose itk and v k to be 
real. 

After minimization, the resulting equations we have to 
solve are 
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K = (ck - nY + K , 
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U k V k 



2 2E k 
A k 

2-Bk ' 



(8) 

(9) 
(10) 

(11) 



with the added constraint that the total density N of 
particles is constant 



N = ($|A>|$) 



d 3 fc 



Pk 



(12) 



III. WAVE FUNCTIONS AND EXCITATION 
SPECTRA 

Numerical solutions are found for an attractive short- 
range interaction, described by a Gaussian potential 
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FIG. 1: (Color online) Scattering length ao and effective 
range of the interaction ro as a function of the potential depth 
Vq/Eq, where the energy scale Eq = fi 2 /(m(r) 2 ). 



V{r) = Vbexp(-6r 2 ). We explore the BCS-BEC 
crossover by fixing the width of the potential (r) = 
2/\^nb and varying the potential depth Vo, or implicitly 
the scattering length ao of the interaction. Fig. ^ de- 
picts all the essential two-body physics of the problem. 
As the depth of the potential is increased from zero, the 
scattering length grows, diverging to negative infinity at 
the first appearance of a two-body bound state. On the 
bound-state side of the Feshbach resonance the scattering 
length re-emerges from positive infinity. Note that this 
procedure does not keep the effective range ro of the po- 
tential fixed. In the low-energy limit, the effective range 
r , like the scattering length ao, can be determined from 
the phase shift 



k cot 5q 



1 1,2 
— + 7j r k 2 , 
ao I 



(13) 



and is shown in Fig. ^ The variation of the effective 
range is small over the region of interest. 

From now on, we shall focus on finite densities and 
the crossover to BCS that can occur in this regime. Like 
previous theoretical studies^, we find that the momen- 
tum distribution pk at constant density smoothly evolves 
from the wave function of a molecule in momentum space 
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FIG. 2: (Color online) Momentum distributions pt (top) 
and the quasi-particle excitation spectra Ek/eF (bottom) as 
a function of k = |k|, at fixed density (&f) and various poten- 
tial scattering lengths do, where ef is the Fermi energy and 
k F (r) « 0.37. 



at large positive (aofep) -1 to a Fermi- like distribution 
at large negative (aokp) -1 , as depicted in Fig. [21 (top). 
During this transition, we see from Fig. [2] (bottom) 
that the quasi-particle excitation spectrum Ek develops a 
minimum at finite momentum, signifying the BCS-BEC 
crossover point defined previously. Because the density 
is non-zero, though small, the crossover point occurs for 
positive aokp, where a two-body molecular bound state 
has already formed. 

The appearance of the minimum in the quasi-particle 
excitation spectrum at non-zero momentum is clearly ev- 
ident in the quasi-particle density of states. As a function 
of momentum, the density of states is defined as 



N k = 



1 k 2 

^2 I dE k 
I dk 



(14) 



Thus, Nk is singular whenever there is a stationary point 
in the excitation spectrum (dEk/dk = 0), except when 
k = 0, where Nk vanishes provided dEk/dk — > no faster 
than k 2 . Fig. [3| shows how Nk varies with scattering 
length at constant particle density. The energy minimum 



for k > shows up as a spike that migrates away from 
k = towards k = kp as (aokp)^ 1 — > — oo, but the spike 
disappears when the energy minimum occurs at k = 0. 
Thus, the BCS-BEC crossover is the point at which a 
spike first appears as we decrease (aokp) -1 . 

A more straightforwardly measurable quantity may be 
the density of states as a function of energy 
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(15) 
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which is depicted in Fig. 0] Approximating the gap pa- 
rameter Afc as a constant, the BCS regime has a singu- 
larity at finite momentum (k = y / 2mp:/h) of the form 



where the amplitude is 
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We see that the size of the singularity will increase the 
further we enter the BCS phase. 

In the BEC limit, the density of states is zero at the 
energy minimum occurring at k = 0, and deep within 
the BEC phase it has the following form close to the 
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FIG. 3: (Color online) Density of states /Vfc (as a function 
of momentum), at fixed fcj? and various potential scattering 
lengths ao, where fcj?(r) « 0.37. 

As well as considering the BCS-BEC crossover at con- 
stant density arranged by tuning the interaction, we also 
examine the density-driven crossover at fixed interac- 
tion. Fig. [S] illustrates the density of states as a func- 
tion of momentum for various particle densities, where 
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FIG. 4: (Color online) Density of states Ne (as a function 
of energy), at fixed fcj?, where kp{r) ~ 0.37, and various 
potential scattering lengths ao. 



FIG. 6: (Color online) Momentum distribution of the conden- 
sate wave function, at fixed density (ftp) and various potential 
scattering lengths ao, where kp(r) ~ 0.37 . 
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FIG. 5: (Color online) Density of states Nk (as a func- 
tion of momentum), at fixed potential scattering length 
{ao/ (r) = 2.4), and various values of fc_F- 



we have chosen ao to be positive so that a molecular 
bound state exists. We observe the same disappear- 
ance of the spike in Nk as {aokw) — > +oo but, unlike 
Fig. El the extreme BCS limit occurs when (aofcp) -1 — > 
since this corresponds to infinite density. Note also that 
both Fig. E2 and [H| show Nk/k — ► const, as (aofcp) 1 — > 
+oo, so that the pair-breaking excitation spectrum has 
(Ek — 5E) (x k 2 , where SE is the binding energy of the 
molecule. This is exactly what we would expect for a 
non-interacting gas of Bose molecules. 

Finally, Fig. shows the behavior of the condensate 
wave function throughout the crossover. Here we see the 
most pronounced features of the evolution from a molec- 
ular state to the BCS limit where pairing exists only on 
the Fermi surface. 
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FIG. 7: Profiles of the chemical potential [i/sf as a function 
of 77 = (aofti?) . The points mark the BEC-BCS crossover. 



IV. PHASE DIAGRAM AND COMPARISON 
WITH ZERO-RANGE POTENTIAL 

A more complete understanding of the BCS-BEC 
crossover can be gained from considering the phase di- 
agram for our finite-range potential. Figure [7] shows 
the variation of chemical potential h/ef as a function of 
(aofc^) -1 , for a set of different scattering lengths ao/(r). 
On each curve we mark the position of the BCS-BEC 
crossover point, which we refer to as the "critical" pair 
(/i c , kp c ). The crossover point always occurs for ao > 
and there is no density-driven crossover when ao < 0. 
Since the size of ao dictates the size of the two-body 
molecular bound state while ao > 0, the crossover point 
in Fig. [7| is driven to lower densities when ao is increased. 
As the zero-density limit (fc^ c =0) is approached, the crit- 
ical chemical potential fx c /£F ~ * 0, which is equivalent to 
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FIG. 8: (Color online) Comparison with Leggett's "univer- 
sal" curve: We keep fc_F constant, and vary the scattering 
length ao- 



FIG. 9: (Color online) Comparison with Leggett's "univer- 
sal" curve: The scattering length ao is held constant, the 
density (kp) is varied. 



the BCS-BEC crossover for a single-parameter potential 
initially discussed by Leggett. However, as fcp c increases, 
the crossover point shifts to negative values of /i c and 
larger (oo^f)" 1 - 

It is useful to compare the results of our calculations 
with previous studies of the BCS-BEC crossover that in- 
volve single-parameter potentials. In the dilute or low- 
energy limit, where kp{r) , {r) / a® <C 1, the effective in- 
teraction U can be written as 



U 



(19) 



and it is independent of the details of the real interac- 
tion. Substituting U for the potential and eliminating 
high energies, Eqs. (|HJ) and l|12|) become, respectively 
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tt r? , (20) 
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where we have introduced the notations fl = ^/sf, A = 
A/sp, and r\ — (agkp) -1 . The above equations, first 
studied by Leggett, are then solved for fi(rj) and A (77). 
For obvious reasons, we will refer to these functions as 
"model independent" or "universal." 

We first compare our results for /j(t7) against Leggett's 
predictions for the case where the scattering length ao is 
varied and the density is held fixed, as shown in Fig|SJ On 
the Leggett curve, the BCS-BEC crossover point (/i = 0) 
occurs at (aokp)^ 1 « 0.55, and our numerical jl(r]) con- 
verges to this universal curve as kp(r) — ► 0. Since typical 
experimental parameters correspond to kp(r) ~ 0.04, the 
ultracold atomic gases clearly lie on the Leggett curve 
and arc well described by a single-parameter potential. 



The value of n/ep for r\ = is usually referred to as 
the "unitarity" Iimi1j2iifliii. In this limit one expects that 
all sensitivity to the detailed features of the interaction 
is lostii^, and the system energy is determined entirely 
by the density. As a result, this limit is particularly sen- 
sitive to many-body correlation effects, and the Leggett 
curve predicts this value as being 0.59 (see alsoii). Re- 
cently, accurate calculations based on the Green's Func- 
tion Monte Carlo (GFMC) methocUi, have lowered the 
upper-bound on this result to (0.44 ±0.01), which shows 
that beyond mean-field effects account for at least a 25% 
improvement in the binding energy over the mean-field 
result. 

Fig. compares the numerical pt(r]) with the universal 
curve for the case where the density is varied and ao/ {r) 
is fixed. Not surprisingly, the numerical p,(r]) approaches 
the universal curve as scattering length increases, and 
we require ao/(r) > 18 before fl(r]) touches the univer- 
sal curve. Furthermore, we see that no matter how large 
we make ao, the universal curve breaks down in the re- 
gion sufficiently close to (aokp)^ 1 = 0. This point corre- 
sponds to infinite density, the BCS limit in the density- 
driven crossover, so we must always have /5(0) = 1, which 
is consistent with our numerical results. If we consider 
p,(rj) for which ao/(r) ~ 72 wc observe that it collapses 
onto the universal curve when {a^kp) -1 > 0.3. This 
translates into the condition kp(r) < 0.046 for the uni- 
versal equations to be valid, which is consistent with 
Fig. El 



V. CONCLUSION 

In this paper we have performed a study of the BCS- 
BEC crossover, at zero temperature, in the single-channel 
model, using a realistic effective interaction. In the mean- 
field approximation the crossover is described using a 
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variational wave function, which interpolates smoothly 
between the BEC and BCS limits. We find that the 
crossover point always lies on the side of the Fcshbach res- 
onance where a molecular bound state exists at zero den- 
sity, and that it converges to that of a single-parameter 
potential in the low density limit. Moreover, our results 
reproduce the usual Leggett picture in the limit of a di- 
lute Fermi gas, as expected. By using a potential of fi- 
nite range, we have shown that the crossover for typi- 
cal experimental parameters in ultracold atomic Fermi 
gases can be described by a zero-range, single-parameter 
potential. We emphasize that a clear signature of the 
crossover point appears in the density of states, which 
could be measured spectroscopically 1 ^. 

The mean-field description of the single-channel model 
presented here is not entirely satisfactory. As described 
above, recent GFMC results provide an upper value for 
the energy density in the dilute limit (77 = 0), which 
represents a 25% improvement over the Leggett limit. 
Therefore, additional work is required in order to include 
next-to-lcading order effects in a quantitative way, and 
study the changes (or lack thereof) in the features of the 
BEC-BCS crossover. Approximations schemes based on 
the two-particle irreducible effective action 1 ^ are readily 
available in quantum field theory where they have been 
recently employed to describe the dynamics of a system 
out of equilibrium, and the dynamics of phase transi- 
tions^. Work is currently under way in order to apply 
such methods to the study of the BEC-BCS crossover 



based on the Hamiltonian QJ. 

The single-channel model is of course a simplification 
of detailed models for the Feshbach resonance which in- 
clude the resonance explicitly, either as a boson (the 
Fermi-Bose modelii) or more generally as three-level — 
or four-level fermi systems. The Fermi-Bose model, as 
well as models with even numbers of fcrmions, all have 
large regimes of parameter space where the resonance is 
sufficiently detuned from the open channel that high en- 
ergy degrees of freedom can be integrated out in favor 
of an effective interaction between two species, recover- 
ing the model discussed in this paper. Close enough to 
resonance the detailed fcrmionic structure of the reso- 
nance may become important. Models where the Fesh- 
bach bound state shares a state with the open channel 
(relevant to 40 K) arc different, and do not reduce to an 
effective single channel model in a simple waj*A&. In fu- 
ture work, we will compare and contrast these different 
models 1 — using a generalization of the variational scheme 
given in this paper. 
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